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Abstract 

In this paper, a de Casteljau algorithm to compute (p, g)-Bernstein Bezier curves based on (p, q )- 
integers is introduced. We study the nature of degree elevation and degree reduction for (p, q)-Bezier 
Bernstein functions. The new curves have some properties similar to (/-Bezier curves. Moreover, we 
construct the corresponding tensor product surfaces over the rectangular domain (u, v) G [0,1] x [0,1] 
depending on four parameters. We also study the de Casteljau algorithm and degree evaluation 
properties of the surfaces for these generalization over the rectangular domain. Furthermore, some 
fundamental properties for (p, (/)-Bernstein Bezier curves are discussed. We get g-Bezier curves and 
surfaces for (u, v) G [0,1] x [0,1] when we set the parameter p 1 = p 2 = 1. 
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1 Introduction 

Recently, Mursaleen et al ® m applied (p, q)-calculus in approximation theory and introduced 
the first (p, g)-analogue of Bernstein operators based on (p, (/[-integers. Motivated by the work of 
Mursaleen et al jjjjj 10], the idea of (p. (/[-calculus and its importance. We construct (p, c/[-Bezier 
curves and surfaces based on (p, (/[-integers which is further generalization of (/-Bezier curves and 
surfaces. For similar works based on (p, (/[-integers, one can refer 1J P2 TTi EE5* H21 [22] • 

It was S.N. Bernstein |Tj in 1912, who first introduced his famous operators B n : CfO, 1] —» C[0,1] 
defined for any n £ N and for any function / G C[0,1] 

Sn(/;*) = E( * )* fc (i-*) n - fc /Q), ^ e [0,1]. (1.1) 

and named it Bernstein polynomials to prove the Weierstrass theorem [?]. Later it was found that 
Bernstein polynomials possess many remarkable properties and has various applications in areas such 
as approximation theory [7], numerical analysis, computer-aided geometric design, and solutions of 
differential equations due to its fine properties of approximation [19] . 


In computer aided geometric design (CAGD), Bernstein polynomials and its variants are used 
in order to preserve the shape of the curves or surfaces. One of the most important curve in CAGD 
[241 is the classical Bezier curve j2] constructed with the help of Bernstein basis functions. 

In recent years, generalization of the Be'zier curve with shape parameters has received continuous 
attention. Several authors were concerned with the problem of changing the shape of curves and 
surfaces, while keeping the control polygon unchanged and thus they generalized the Bezier curves in 

mmm- 


The rapid development of g-calculus [23] has led to the discovery of new generalizations of Bern¬ 
stein polynomials involving g-integers [HUHES] • The aim of these generalizations is to provide 
appropriate and powerful tools to application areas such as numerical analysis, computer-aided geo¬ 
metric design, and solutions of differential equations. 

In 1987, Lupa§ [8] introduced the first ^-analogue of Bernstein operator as follows 


Ln,q(f',x) — ^ ' 
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fc(fc—i) 


q 2 x k (1 — x) 
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and investigated its approximating and shape-preserving properties. 


( 1 . 2 ) 


In 1996, Phillips ]T7] proposed another q -variant of the classical Bernstein operator, the so-called 
Phillips g-Bernstein operator which attracted lots of investigations. 


B n , q (f;x) = ^2 
k—0 
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n—k —1 

x k n (1 -q s x) 

9 s=0 
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where B n ^ q : C[0,1] —>• C[0,1] defined for any n G N and any function / G C[0,1]. 


(1.3) 


The q -variants of Bernstein polynomials provide one shape parameter for constructing free-form 
curves and surfaces, Phillips g-Bernstein operator was applied well in this area. 

In 2003, Oruk and Phillips [HI used the basis functions of Phillips g-Bernstein operator for 
construction of g-Bezier curves, which they call Phillips g-Bezier curves, and studied the properties of 
degree reduction and elevation. 

Thus with the development of (p, g)-analogue of Bernstein operators and its variants, one 
natural question arises, how it can be used in order to preserve the shape of the curves or surfaces. 
In this way, it opens a new research direction which requires further investigations. 


Before proceeding further, let us recall certain notations of (p, g)-calculus . 


The (p, g) integers [n] Pj9 are defined by 

v n — q n 

[n]p,q= -, n = 0,1,2, p > q > 0. 

p-q 

The formula for (p, g)-binomial expansion is as follow: 
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where (p, g)-binomial coefficients are defined by 


n 
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p,q [^’]p,9-[ n ^]p,9' 

Details on (p, g)-calculus can be found in (6j [5J [21]. 

The (p, q)-Bernstein Operators introduced by Mursaleen et al is as follow: 
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Note when p = 1, (p, g)-Bernstein Operators given by (11.41) turns out to be q-Bernstein Operators. 
Also, we have 
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Again by some simple calculations and using the property of (p, g)-integers, we get (p, g)-analogue 
of Pascal’s relation as follow: 
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(1.5) 

( 1 . 6 ) 


We apply (p, g)-calculus and introduce first the (p, (?)-Bezier curves and surfaces based on 
(p, ^-integers which is further generalization of g-Bezier curves and surfaces, for example, [181 iT5]. 

The outline of this paper is as follow: Section 2 introduces a (p, g)-analogue of the Bernstein 
functions B and their Properties. Section 3 introduces degree elevation and degree reduction 
properties for (p, g)-analogue of the Bernstein functions. Section 4 introduces a de Casteljau type 
algorithm for B k '™. In Section 5 we define a tensor product patch based on algorithm 1 and its 
geometric properties as well as a degree elevation technique are investigated. Furthermore tensor 
product of (p, g)-Bezier surfaces on [0,1] x [0,1] for (p, q)-Bernstein polynomials are introduced and 
its properties that is inherited from the univariate case are being discussed. 


2 (p, g)-Bernstein functions 

The (p, (?)-Bernstein functions is as follows 


s k P :(t) = 
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( 2 . 1 ) 


3 
































where 


n—k—1 

(i-t);~ k = n (p s -9 s t) 

s =0 

2.1 Properties of the (p, g)-analogue of the Bernstein functions 

Theorem 2.1 The ( p , q)-analogue of the Bernstein functions possess the following properties: 

(1.) Non-negativity: Bp'™(t) > 0 k = 0,1, ...,n, t £ [0,1], 

(2.) Partition of unity: 

n 

b p,q W = 1 ’ f° r ever y t e [°> 

k =0 


(3.)Both sided end-point property: 

B k p %( o) = 
i& n (i) = 


if k = 0 

if k = n 
k ^ n 


when p = 1, i/ien both side end point interpolation property holds. 

(3.) Reducibility: when p = 1, formula (2.1) reduces to the q-Bernstein bases. 


Proof: All these property can be deduced easily from equation m- 
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q = 0.8 and p = 1 


Figure 1: ‘Cubic Bezier blending functions’ 



q = 0.5 and p = .8 


Figure 2: ‘Cubic Bezier blending functions’ 
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Fig. [2]shows the (p, g)-analogues of the Bernstein basis functions of degree 3 with q = 0.5, p = .8. 
Here we can observe that sum of blending fuctions is always unity and also end point interpo¬ 
lation property holds, when we put p = 1, it turns out to be g-Bernstein basis which is shown in Fig. [lj 

Apart from the basic properties above, the (p, g)-analogue of the Bernstein functions also satisfy 
some recurrence relations, as for the classical Bernstein basis. 


3 Degree elevation and reduction for (p, q)- 
Bernstein functions 

Technique of degree elevation has been used to increase the flexibility of a given curve. A degree 
elevation algorithm calculates a new set of control points by choosing a convex combination of the 
old set of control points which retains the old end points. For this purpose, the identities (13.31) . (13.dl) 
and Theorem (13.21) are useful. 


Theorem 3.1 Each (p, q)-Bernstein functions of degree n is a linear combination of two {p,q)- 
Bernstein functions of degree n — 1 : 
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Proof: On using Pascal’s type relation based on (p, g)-integers i.e (11.51) . we get 
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Degree elevation 
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By some simple calculation, we have 
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finally we get 
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Theorem 3.2 Each (p, q)-Bernstein function of degree n is a linear combination of two ( p,q)- 
Bemstein functions of degree n + 1. 
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Proof From equation 13.31 3.41 we can easily get 
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For n = 3, the blending functions are given by: 
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We observed that the both sided end point interpolation property and partition of unity property 
always holds in case of (p, g)-Bernstein functions. 


The de Casteljau algorithm describes how to subdivide a Bezier curve, when a Bezier curve is 
repeatedly subdivided, the collection of control polygons converge to the curve. Thus, the way of 
computing a Bezier curve is to simply subdivide it an appropriate number of times and compute the 
control polygons. 


4 (p, t/)-Bernstein Bezier curves: 

Let us define the (p, g)-Bezier curves of degree n using the (p, <?)-analogues of the Bernstein functions 
as follows: 


P(*;p, «) = £?**&(*) (4-i) 

i=0 

where Pi G R 3 (i = 0, l,...,n) and p > q > 0. Pi are control points. Joining up adjacent 
points Pi , i = 0,1,2,..., n to obtain a polygon which is called the control polygon of (p, g)-Bzier curves. 


4.1 Some basic properties of (p,q )~Bezier curves. 

Theorem 4.1 From the definition, we can derive some basic properties of (p, q)-Bezier curves: 

1. (p, q)-Bezier curves have geometric and affine invariance. 

2. ( p,q)-Bezier curves lie inside the convex hull of its control polygon. 

3. The end-point interpolation property: P(0;p, q) = Po, P(l;p,< 7 ) = P n - 
3. Reducibility: when p = 1, formula \4.1\ gives the q-Bezier curves. 

Proof. These properties of (p, g)-Bezier curves can be easily deduced from corresponding prop¬ 
erties of the (p, q)-analogue of the Bernstein functions. 
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4.2 Degree elevation for (p, q)- Bezier curves 

(p, q )-Bezier curves have a degree elevation algorithm that is similar to that possessed by the classical 
Bezier curves. Using the technique of degree elevation, we can increase the flexibility of a given curve. 
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The statement above can be derived using the identities (l3.3l) and (13.4H . Consider 
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where Pljis defined as the zero vector. Comparing coefficients on both side, we have 
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where r = 1,.n and i = 0,1 , 2, ..., n + r 

When p = 1, formula 14.41 reduce to the degree evaluation formula of the g-Bezier curves. If we 
let P = (Po, Pi,..., P n ) T denote the vector of control points of the initial (p, g)-Bezier curve of degree 
n,andP( 1 ) = (Po*,P 1 * 1 ... 1 P* +1 ) the vector of control points of the degree elevated (p, g)-Bezier curve 
of degree n+1, then we can represent the degree elevation procedure as: 


P (1) = T„ +1 P 


where 
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For any l € N, the vector of control points of the degree elevated (p, g)-Bzier curves of degree n + Z 

is: PO) = T n+ i T n+2 . ...T n+ iP. As l —>• 00 , the control polygon P (1 ' converges to a (p, g)-Bezier 

curve. 


4.3 de Casteljau algorithm: 

(p, g)-Bezier curves of degree n can be written as two kinds of linear combination of two (p, g)-Bezier 
curves of degree n — 1, and we can get the two selectable algorithms to evaluate (p, g)-Bezier curves. 
The algorithms can be expressed as: 

Algorithm 1. 


P?(t;p,q) = P? = Pi * = 0,1,2. ,n 

P r i (t-,p,q)=p r -H Pl~l(t-,p,q) + (dV^- 1 - q r -H) P^ 1 (t, p, q) 
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r = l,..., n, * = 0,1,2. ,n — r., 

Then 

n —1 

p(*;p,<?) = = B pT r tf) = = p S 

i=0 

It is clear that the results can be obtained from Theorem (??). When p = 1, formula (14.31) and 
(14.41) recover the de Casteljau algorithms of classical g-Bezier curves. Let P° = (P 0 , Pi,..., P n ) T , 
P r = (Pp,PC,....,P^_ r ) T , then de Casteljau algorithm can be expressed as: 

Algorithm 2. 

P r (t;p,q) = M r {t\p,q)....M 2 {t-,p,q)M 1 {t\p,q)P° (4.6) 

where M r (t\p , q) is a (n — r + 1) x (n — r + 2) matrix and 
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5 Tensor product (p, f/)-Bernstein Bezier surfaces on 

[0,1] x [0,1] 

We define a two-parameter family P(u, v) of tensor product surfaces of degree m x n as follow: 


p («, *) = E E P ( M ) («). *0 e [0,1] X [0,1], (5.1) 

2=0 J =0 

where P,;.j el 3 (i = 0,1, = 0,1, ...,n) and two real numbers pi > q\ > 0, P 2 > qi > 0, 

hpi'qi ( M ) ; b J p™ q2 ( v ) are (p, g)-analogue of Bernstein functions respectively with the parameter pi,qi 
and P 2 , q 2 • We call the parameter surface tensor product (p, g)-Bezier surface with degree mx n. We 
refer to the P, j as the control points. By joining up adjacent points in the same row or column to 
obtain a net which is called the control net of tensor product (p, g)-Bezier surface. 


5.1 Properties 

1. Geometric invariance and affine invariance property: Since 

m n 

EE«» = 1. (5-2) 

2=0 j =0 

P(w, v) is an affine combination of its control points. 

2. Convex hull property: P (u,v) is a convex combination of P t j and lies in the convex hull 
of its control net. 

3. Isoparametric curves property: The isoparametric curves v = v* and u = u* of a tensor 
product (p, g)-Bezier surface are respectively the (p, g)-Bezier curves of degree m and degree n, namely, 


m 

■*) = E ( 

f n \ 

E p m^>*) 


u e [0,1]; 

i=0 

^ 3=0 J 



II 


I 

, v e [o, l] 

3=0 

' z=0 * 




The boundary curves of P(u,u) are evaluated by P(xi, 0), P(u, 1), P(0,u) and P(l,v). 

4. Corner point interpolation property: The corner control net coincide with the four 
corners of the surface. Namely, P(0,0) = Po.cn P(0, 1) = Po,m P(l> 0) = P mi o, P(l, 1) = P m ,n, 

5. Reducibility: When pi = P 2 = 1, formula (EED reduces to a tensor product g-Bezier patch. 
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5.2 Degree elevation and de Casteljau algorithm 

Let P(u,v) be a tensor product {p,q )~Bezier surface of degree m x n. As an example, let us take 
obtaining the same surface as a surface of degree (m + 1) x (n+ 1). Hence we need to find new control 
points P* • such that 


m n m+1 n+1 

p(v)=EE p y = EE p, i. j b£+Hv) ( 5 . 3 ) 

2=0 j= 0 2=0 j — 0 


Let a.i=p 1 ‘ 
Then 
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(5.4) 


which can be written in matrix form as 
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The de Casteljau algorithms are also easily extended to evaluate points on a (p, q )-Bezier surface. 
Given the control net Pjj G R 3 , * = 0,1,..., to, j = 0,1,..., n. 


P"j>,*0 = P ij = p i,j * = 0,1,2., m; j=0,l,2...n. 


p- j r («,*o = [bbr^ 1 - 

r = l,...,fe, fc = min(m, n) * = 0,1,2., to — r; j = 0,1, ....n — r 


1 

f- -3 

1 

1 

pr— l,r—1 

^i,j +1 , 


Y( j r—j —1 r —1 

( q 2 p 2 v ) 

pr-l,r-l 
. *+l,J 

Tjf— 1 ,r—1 
^i+l,j+l _ 


r— 1 

L P 2 v J 


(5.5) 


or 


p -j°(«, *0 = p °j° = p ij * = 0,1,2., to; j = 0,1, 2.,.n. 


-1 

1 

1 

pr— l,r—1 

^i,j+ 1 


["/ r— 1 j r—j —1 \1 

\p 2 -P2Q2 v ) 

pr-l,r-l 

. i+l,j 

13 r — 1, r —1 
^i+i,j+i _ 


7 r—j —1 

l P2Q2 v J 


pff(u,«) = [br 1 - Pbr i_ b) 

r = l,...,fc, fc = min(TO,n) * = 0,1,2., to — r; j = 0,1, ....n — r 

(5-6) 

When to = n, one can directly use the algorithms above to get a point on the surface. When to n, 
to get a point on the surface after k applications of formula (15.51) or (15.61) . we perform formula (14.61) 
for the intermediate point Pb . 

Note: We get g-Bezier curves and surfaces for (u,v) G [0,1] x [0,1] when we set the parameter 
Pi = P 2 = 1 as proved in mm- 
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6 Shape control of (p, g)-Bernstein curves 

We have constructed ( p , g)-Bernstein functions which holds both the end point interpolation property 
as shown in figure Q] and [2j Parameter p and q has been used to control the shape of curves and 
surfaces: ifO<g<p<l, asp and q decreases, the curve moves close to the control polygon, as p 
and q increases, the curve moves far away from the control polygon; If p > 1 and q > 1, the effects of 
p and q are opposite, as p and q decreases, the curve moves far away from the control polygon, as p 
and q increases, the curve moves close to the control polygon. 

Figure 1314151 Gl mdTl shows (p, g)-Bernstein function approximating the surface generated by control 
points. 



Figure 3: ‘Bezier curve 



Figure 4: ‘Bezier curve 
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Figure 5: ‘Bezier curve’ 



Figure 6: ‘Bezier curve’ 
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